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SUMMARY 


A numerical study of the effects of viscous— inviscid inter- 
action in three-dimensional duct flows is presented- In particu- 
lar we consider interacting flows for which the oncoming flow is 
not fully-developed- In this case there is a thin boundary layer 
still present upstream of the of the surface distortion, as 
opposed to the fully-developed pipe flow situation wherein the 
flow is viscous across the entire cross section. 
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1. INTRODUCTION 


There exist several reliable methods (e.g. Vatsa & Davis 
1973, Blottner 1975, Cebeci , Kaups & Ramsey 1977, McLean & 

Randall 1979, Lindhout, Moek, Deboer & van der Berg 1979) -for 
calculating standard (or direct) three-dimensional boundary— 1 ayer 
flows with the pressure prescribed- However, such methods are 
seldom useful in real flows because of the classical separation 
singularity, among other things- Instead vi scous-i nvi sci d inter- 
action generally occurs with both the pressure and the boundary- 
layer displacement unknown in advance, just as in two-dimensional 
f 1 ows. 

Smith (1983) has studied the properties of, and a finite- 
difference approach for, interactive three-dimensional boundary 
layers in external flow problems. He found that, contrary to the 
results of two-dimensional interactions, the three-dimensional 
inverse problem is an elliptic one and therefore standard foward— 
marching in the flow direction, e.g. through cross-flow planes, 
is strictly i nappropr i ate/dangerous because an explosive 3-D free 
interaction can be initiated. The same phenomena can take place 
in the interactive case, or whenever the pressure is not de- 
scribed. To overcome these problems. Smith (1983) and Edwards, 
Carter, & Smith (1987) developed an unconventional finite- 
difference approach for the external three-dimensional inverse 
and triple-deck problems which addresses the major numerical 
difficulties of ellipticity and violent three-dimensional free 
interactions which are inherent in the standard f oward-marchi ng 
schemes. In this paper the numerical scheme developed by Smith 


1 



for three-dimensional interactive flows is generalized to include 
three-dimensional internal flows with vi scous— i nvi sci d interac- 
tion, such as in a transition duct- More specifically, we study 
the interactive flow field present in a duct when the oncoming 
flow is not yet f ul ly— developed , i-e-, in those situations where 
the entry length isn't long enough for the flow to become fully- 
developed. In these cases there will still be a "thin" boundary 
layer on the wall of the duct and hence the possibility of 
viscous-invi scid interaction in the vicinity of the surface dis- 
tortion with a nonzero incremental displacement thickness 
function , A (x , €0 « 
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2. GOVERNING EQUATIONS 


For a developing -flow in the geometry shown schematically in 
Figure 1, it can be shown, -following Smith (1976), that the flow 
structure is three— tiered, consisting of a viscous sublayer I of 
thickness 0(Re~ 2/5 ), a main deck II of thickness 0(Re _1/5 >, 
which is an inviscid rotational perturbation of the incoming 
boundary-layer flow, and a potential core region III of thickness 
0(1), where Re = a*U*/V is a representative Reynolds number based 
on the radius a* of a straight pipe of circular cross section, 
the centerline velocity U* in the oncoming flow, and the fluid's 
kinematic viscosity, V. Furthermore, we assume that the distur — 
bance (e.g. corner, hump, indentation, etc) size is of the same 
order of magnitude as the perturbation to the oncoming shear 
flow. Unlike most problems involving vi scous-invi sci d interac- 
tions, the problem considered here reduces to a consideration of 
the equations of motion in the viscous sublayer I and the outer 
core region III. Thus in I we let 


r = 1 - 

Re" 2/5 , Y = 0(1), 

x = od), e = od) 

u(r,©,X) 

= Re“ 1/5 U(Y,0,X) 

^ — — — , 

v(r ,©,X) 

= Re“ 3/5 V(Y,©,X) 


w(r , ©, X ) 

= Re _1/5 W(Y,©,X) 

+ ... , 

p(6,X> = 

Re" 2/5 P(e,X) + . 

* ■ f 


where (r,0,X) are suitably nondimensional ized cylindrical polar 

V V 'V 'V 

coordinates, (v f w,u> the corresponding velocity components and p 
the nondimensional pressure- Then it is readily shown that in 
region I the Navi er-Stokes equations reduce to the three- 
dimensional boundary-layer equations, viz: 
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U x + V Y + W e = O 



(2-2) 


uu x + vu Y + wu e = -p x 

+ Uyy , 


(2-3) 


uw x + vw Y + ww e = -P 0 

+ Wyy , 


(2-4) 

with 

boundary conditions 





U = V = W = 0onY = 

F(X,9) , 


(2-5) 

where F(X,9) is a suitably normalized 

wall shape function 

such 

that 

F — i- 0 as X — -oo, and F(X,€h-2L) 

= F(X,6) ; 



U(Y,0,X) Y, V,W,P, 

, A — ■+ 0, as 

X — -oo. 

(2 

and 






U~Y + A(X,e), Y-> 

oo, all X, 


(2-7) 


W x -P 0 /(Y+A) , Y -s- 

oo, all X , 


(2-8) 


and periodicity in 6. 

Here (2-5) corresponds to the no-slip condition, (2-6) corres- 
ponds to matching with the upstream incident -flow, and (2-7) 
results from a matching of the boundary-layer flow in I to the 
inviscid rotational flow in region II. The function — A(X,0) 
represents the unknown incremental displacement thickness and is 
linked to the unknown pressure P(X,6) by means of the quasi - 
in viscid properties in the core flow III. 

In the core region we must solve for P^ ny (r,©,X), where 

V^inv = 0, 8 (2-9) 

along with the boundary conditions 

P inv = P(X,€» @ r = 1, 

< P inv>r - "AXX ®r=l, (2-10) 

regularity at r = 0. 

Finally, to facilitate the numerical computations the 
governing equations in the viscous sublayer are rewritten using 
Prandtl 's transposition theorem, viz 
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y » Y - F(X,e> , x = X, 

v = V - UF X - WF e , (2-11) 

u(y,0,x) = U(Y,e,X), w(y,0,x) = W(Y,©,X), 


p(x,©> = P(X,©) . 

Then equations (2-2) - (2-4) become 

u x + v y + w e = 0, (2-12) 

UU X + vu y + wu e = -P>; + u yy» (2-13) 

uw x + vw y + WW 6 = -p e + w yyt (2-14) 

with boundary conditions 

u = v = w = 0 on y = 0, all x , 0, (2-15) 

u ~ y + F (x , 8) + A(x,8), y — *- oo, all x, (2—16) 

w x ~ -p e /(y+A+F>, y oo, all x, (2-17) 

u ~ y, v , w , p , A 0 , x — -oo, all y, (2-18) 

and periodicity in O. 


As noted by Smith (1976,1980) these equations are elliptic 
in the pressure field and therefore require a modification of the 
multi— sweep f orward-marching methods developed successfully for 
two-dimensional flow problems. In the present work the ellip— 
ticity embedded in the pressure field needs special attention for 
otherwise an explosive 3D free interaction can be set off in a 
forward-marching procedure. Consider the sum of the x— derivative 
of (2-13) , the y-derivative of (2-14) and the x— derivative of (2— 
12). then defining the skewed shears in the following manner: 

u = u x + w $ , v = v x , (2-19) 

the fundamental governing equations can be written as 


U X + v e = °» 

(2-20) 


uu x + U y v = -E + u yy - S 

(2-21) 


where 
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( 2 - 22 ) 


S(x,y,©> = u 2 + vu sy + <wu e ) x + u e w x 

+(vw y ) e + (ww e ) e , 

and 

E = Pxx + Pee* (2-23) 

In these equations u, v, and E are treated as the unknowns and 
the -following boundary conditions are applied 

u = v = 0@y = 0, all x , (2-24) 

□ , v, E — > 0 , -x — oo, (2-25) 

u - (A+F) xf v+E - u (A+F) xx - (A x +F x ) 2 , y oo, all x. (2-26) 
Note that in these variables the governing equations (2-20) - (2— 
21) are quasi two-dimensional and linear in the unknown 
vari abl es. 

In the core region. III, we must solve 

= E r|r + r~ l E r + r~ 2 E ee + E xx = O, (2-27) 

where 


E(x,r,6> = p xx + p ee 

with 

I(x,©) = E(x,l,©) , 

and 


(2-28) 


(2-29) 


E(x,r,6) bounded as r —*■ 0. (2-30) 

s 

In the numerical computations we -find it convenient to 
introduce an alternative variable to the displacement -function 
A(x,6), namely B(x,0), where 

B (x , 9) = A xx . (2-31) 

The (key) matching relationship is then given by 

E r (x,l,©) = — (B xx + B 0e ). (2-32) 

We also -find it convenient to rewrite (2-26) somewhat, using 
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(2-31) together with (2-16), to give 

x 

G - J Bdx + F x , (2-33) 

— CD 

v + £ ~ -u(B+F xx > - v£ (2-34) 

as y — ■* oo -for all x, ©. 

For convenience, symmetry about the plane © = O was taken 
with reflection conditions (u,v,w,p) (x ,— y ,6) = (u, v ,— w ,p> (x ,y , ©> ; 
to account for periodicity, reflection conditions were also 
applied at © = L , for some constant L. 
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3. NUMERICAL PROCEDURE 


The nature of the boundary-layer problem described above 
suggests an iterative multi— sweep technique, using forward march- 
ing in a quasi two-dimensional manner to solve for u, v and E 
with u, w, and S assumed known from the previous iteration (or an 
initial guess). Simultaneously, we use standard central differ- 
ences in order to solve (2-27). 

Specifically, given a guess or update for u,v,w and there- 
fore S everywhere, (2—20) - (2-21), (2—24) - (2-26), are marched 
forward in x, while at the same time (2-27), (2-29), (2—30) are 

solved along a line of varying r; this then determines u,v,E,E,B 
(and. hence A) along a line of constant spanwise coordinate. The 
process is then repeated at all other spanwise coordinates. Next 
the surface pressure p(x,0) is found from the Poisson equation 
(2-33). Updated values of w are found next by marching (2-14) 
forward in x for each ©, subject to (2—15) - (2—17) , assuming 
that u,v,p and w^ are known. Finally updated values of u and v 
are found by an integration of (2-19) with respect to x along 
with (2-15) and (2-16). All the above steps constitute one 

global iteration. Convergence is attained when a global convei 

s 

gence test on u is satisfied. 

The main features of the numerical scheme are the following. 
Two and three-point differencing in y is used for (2-20) and (2- 
21), respectively, with (2-26) applied at y = y^, two-point 
differencing in x and three-point differencing in © when S is 
evaluated, at a given x,© location. This scheme reflects the 
fact that the flow is generally parabolic in x but elliptic in y. 
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Three— point, central differencing is used to approximate (2—27) in 
all three dimensions , and (2-32) is approximated by three-point 
one-sided differencing. Supposing we have n points in y, m points 
in r, then at each x,0 station the difference approximation of 
(2—20), (2—21) and (2-27), together with the interface conditions 

(2—32) — (2—33), may be written schematically in the following 
form (where an X denotes a non-zero entry and 0 a zero entry) : 


•e 


* 



xxoo 

0 

(2-20) 


v 2 

xxxx 

X 

(2-21) 


u 2 

xxxxoo 

0 

(2-20) 


v 3 

xxxxxx 

X 

(2-21) 


u T 

xxxxoo 

0 

(2-20) 


v 4 

xxxxxx 

X 

(2-21) 


“4 

v 5 

xxxxoo 

0 

(2-20) 


xxxxxx 

X 

(2-21) 


,. r u 

u 5 

* 

(2-20) 

(2-21) 

(2-20) 

(2-34) 

(2-33) 

(2-32) 

(2-27) 

m m 

XXXXOOOO 

XXXXXXOX 

xxxxoo 

xoxx 

xxo 

xxxx 

XXX 

i 

0 

m 

m 

In 

u n 

B 

?=E 1 

E 2 

E 3 

E 4 

(2-27) 

(2-27) 

(2-27) 

XXX 

• ■ • ■ 
m 

xxxo 

XX 


• 

m 

k-2 

E m- 1 

■ 


m 




Here equation numbers denote* the equation which the row 
approximates. Notice that 

= Gj = 0 , E b = 0 . (3-2) 


The solitary column of entries in the upper part of the 
matrix corresponds to the E term in (2—21). (3—1) is then amen- 

able to standard Gaussian elimination procedures. (2-33), which 
determines the surface pressure, was also approximated by central 
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differences, and was solved by Gauss-Seidel sweeps (without 
pivoting) , with an over-relaxation parameter until a prescribed 
tolerance is met. The solution for w is similar to that for u 
and v except that only a tri-diagonal matrix requires inversion 
(performed again using Gaussian elimination). Finally the x- 
integrations to compute u and v were performed using the trape- 
zoidal rule. The overall scheme is nominally second-order accur 
ate in the mesh widths Ax, Ay, A©, and Ar. 

A crucial feature of the scheme developed as parrt of this 
research effort is the inherent coupling between the viscous 
boundary-1 ayer solution and the inviscid core solution which is 
carried out simultaneously in the spirit of the scheme proposed 
by Veldman (1979) in treating two-dimensional incompressible 
flows, using the Hilbert integral approach. This approach has 
also recently been used by Bodonyi & Duck (19S7) to succesfully 
treat three-dimensional external interacting flows. 
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4. RESULTS 


The numerical method discussed in the previous section has 
been applied to several test problems. The model problem 
corresponds to that of a nonf ul ly— developed flow upstream of the 
surface distortion given by a bounded transition duct with F(x,6) 
defined by 

hx 3 ( 1 + x 3 ) ~ 1 cos (26) , x > 0, 

F(x,e) = (4-1) 

0. x < 0. 

and h ranging between 0.6 and 3.2. 

In this case the duct undergoes a transition from that of a 
circular cross section duct far upstream to that of an elliptical 
cross section far downstream, as shown schematically in Figure 2. 

Solutions for the fully developed flow problem upstream of 
the surface distortion wherein A(x,0) = O have also been found 
for the same duct shape, ince the solutions for both of these 
problems are qualitatively similar, only the results for the 
interactive case are given here. Figures 3 -7 give the axial 
and azimuthal skin friction distributions and the wall pressure 
distribution along the line of symmetry for several values of h. 
No flow reversal was encountered in the axial direction for 
values of h considered in this report, although the azimuthal 
skin friction Wy does exhibit a reversal of sign sufficiently 
far downstream for h > 1. All efforts to extend the computations 
to larger values of h such that axial flow reversal occured were 
unsuccessful. The difficulties in obtaining solutions for h > 

3.2 are apparent in the pressure distribution shown in Figure 7. 

A possible reason for this is that the governing equations are 
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•treated in a quasi -two-dimensional forward marching manner which 
may not properly account -for the transverse velocity component 
when it becomes sufficiently large as in the bounded transition 
duct when h > 3.2. Unfortunately, within the time frame of this 
research effort this difficulty was not resolved. 

A numerical study of noncircular incident ducts has also 
been considered. In this case the incident skin friction depends 
on 8, e.g. parallel flow through an ellipse. The governing 
equations are the same as before except that now the outer 
boundary condition (2—16) is replaced by 

u ^ X(8)f y + A(x,8> + F (x , 8) } . (4-2) 

Here the parameter X(8) gives the azimuthal distribution of the 
incoming shear flow in the viscous sublayer u = X(8)y. 

The numerical scheme discussed above has been adapted to 
include this effect and representative results for the skin 
friction distributions when A(x,8) = 0 and F(x,8)is given by (4- 
1) are shown in Figure B for h = 2 and X(8) — 0.25 + sin 2 8. 

Finally, as part of this research effort similarity 
solutions for x — oo of equations (2-12) — (2—18) for an 
“unbounded" transition body shape of the form 

F(x,©> ~ x 1 /3 cos (20) * (4-3) 

have been investigated with the goal of gaining a better 
understanding of the flow structure far downstream in the duct. 

By expanding in powers of 8, the governing equations can be 
reduced to a set of ordinary differential equations, but in the 
process several undetermined parameters are introduced. It was 
found that the appropriate similarity equations permit multiple 
solutions, as well as regions in parameter space where no 
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solutions could be -found depending on the range of values o-f the 
parameters mentioned above. For -further details concerning these 
similarity equations, the interested reader is referred to 
Bennett (1986) who studied these similarity equations in another 
context. 
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Figure 1. Schematic of the Flow Field Geometry. 


15 



INCOMING 

FLOW 



SECTION 


Figure 2. Schematic of Duct Geometry for Bounded 
Transition Flow Given by Equation 36. 
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Figure 3. Axial Wall Shear Stress for a Bounded Transition Duct 


P(x,t*) 
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Figure 4. Pressure Distribution for a Bounded Transition Duct, h = 1,0 
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Figure 5. Axial Wall Shear Stress 



w y (x,0,£) 
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Figure 6. Azimuthal Wall Shear Stress for a Bounded Transition Duct 


P(X,n>) 
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Figure 7. Pressure Distribution for a Bounded Transition Duct 



